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Abstract 

At the quark-hadron transition, when quarks get confined to hadrons, certain 
orbitally excited states, namely those which have excitation energies above the 
respective L = states of the same order as the transition temperature T c , 
may form easily because of thermal velocities of quarks at the transition tem- 
perature. We propose that the ratio of multiplicities of such excited states 
to the respective L = states can serve as an almost model independent 
signal for the quark-gluon plasma formation in relativistic heavy-ion colli- 
sions. For example, the ratio R* of multiplicities of D^j(2317)(J p = + ) 
and L>^ ± (2112)( J p = 1~) when plotted with respect to the center of mass 
energy of the collision ^fs (or vs. centrality/number of participants), should 
show a jump at the value of yfs beyond which the QGP formation occurs. 
This should happen irrespective of the shape of the overall plot of R* vs. yfs. 
Recent data from RHIC on A* /A vs. N par t for large values of N par t may be 
indicative of such a behavior, though there are large error bars. We give a 
list of several other such candidate hadronic states. 
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I. INTRODUCTION 

There are strong reasons to believe that in ultra-relativistic collisions of heavy nuclei, 
a hot dense region of quark-gluon plasma (QGP) may get created. There is a wealth of 
data which strongly suggests that this transient QGP state may have been achieved in these 
experiments [1] . Still a conclusive evidence for the QGP detection is missing. Among many 
signals which have been proposed for the detection of QGP, probably the suppression of 
J/ijj had the potential of being the cleanest signal [2]. However, due to various uncertainties 
involving final state interactions etc. it has not been possible so far to use even this signal 
to make definitive statements about QGP detection [3]. Most other signals depend more 
sensitively on the details of plasma evolution as well as on the details of the quark-hadron 
transition dynamics, and hence the predictions become crucially dependent on underlying 
models used to describe the dynamics of various stages of the parton system evolution [4]. 

Ideally one would like to have a signal which depends least on the details such as the 
order of quark-hadron transition etc. In view of lattice results, which suggest that quark- 
hadron transition may be a cross-over [5], one would like to only use the single defining 
property of the quark hadron transition, namely that quarks and gluons are deconfined 
above the transition temperature T c , while they are confined into hadronic states below 
T c . (For a smooth cross-over, there may not be a sharp value of T c at which confinement 
properties of quarks change abruptly. In such a case, by the transition temperature T c we 
will mean a range of temperature near the value T = T c across which the confinement of 
quarks significantly changes. If this range is narrow enough then use of a definite value 
for T c will still capture the essential physics of our model.) In the following we argue 
that the ratio R* of multiplicities of certain specific orbitally excited hadronic states to the 
respective L = states, where L is the orbital angular momentum, may provide such a 
model independent signal for QGP formation in relativistic heavy-ion collisions. (In the 
following, by orbitally excited state, and orbital excitation energy, we will always mean the 
L = 1 orbital excitation above the corresponding L = state, and the energy splitting of 
these two levels.) For example, the ratio R* of multiplicities of DJj(2317)(J p = + ) and 
D* S ± (2112)(J P = l - ) when plotted with respect to the center of mass energy of the collision 
V^s, or centrality /number of participants (N pari ), should show a jump at the value ,Js~ c of 
yt~s (or centrality) beyond which the QGP formation occurs. 

It is important to clarify here that we do not discuss how R* should behave as a function 
of y/s (or, N part ). That is a separate issue and involves factors like the final state interaction 
etc. Our only claim is that, whatever be the form of the curve of R* vs. ^J~s (increasing or 
decreasing), there should be a jump enhancement at (or near the value) ^fsl (or, at some 
specific value of N part ) beyond which QGP formation happens. In this context we note 
that the recent data from RHIC on A*/A shows a strong decrease as a function of N part 
[6,7]. However, as just mentioned above, this is of no consequence to our model. What is 
important is to look for any possible jump in this decreasing curve at some specific value of 
N part . We note from ref. [6] that the value of the multiplicity ratio for A* /A for very large 
Npart (beyond 350) may be consistent with a slight increase, though completely within the 
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error bar, to conclude anything. (This possible increase is less clear in ref. [7], which could 
be due to different binning.) What is important is that if there is any rapid increase in 
R* beyond some N part then it will be very hard to explain it with existing models where a 
decrease in the multiplicities of resonances arising from final state interaction in the hadronic 
medium is naturally expected. (Note that any resonance regeneration should show smooth 
dependence on y/s.) As we will see below, our model will naturally account for any such 
increase. 

We give a list of several hadronic states which could provide such a signal for QGP 
formation in relativistic heavy-ion collisions. We mention that in order that these enhanced 
multiplicities of orbitally excited states survive up to the hadronic chemical freezeout, the 
chemical equilibrium should not last for long after the hadronization. This may not be 
an unreasonable assumption as chemical freezeout temperatures are close to the expected 
value of the transition temperature [8] . It is also expected that resonances may not achieve 
chemical equilibrium. 

The paper is organized in the following manner. In section II, we describe the physical 
picture underlying the signal we are proposing. In section III, we discuss the details of the 
signal and certain special features specific to this signal. In section IV we discuss a rough 
estimate of the ratio R* based on a recombination model. In section V we give a list of 
candidate hadronic states which can be used for this signal. Section VI gives conclusions. 

II. PHYSICAL PICTURE OF THE MODEL 

We start by trying to form a physical picture of the process of hadron formation at the 
transition stage. We will restrict our attention to hadrons with not very high Pp, so that 
the hadronization will be expected to be dominated by something like the recombination 
process, rather than by fragmentation [9,10]. This may cover hadrons starting from very 
low P T , up to hadrons with P T as large as few GeV/c, see ref. [9]. (As we will later argue, 
there are other arguments in our model suggesting that the signal we propose will be more 
prominent when restricted to low Pp hadrons.) Further, as in certain recombination models 
[9], at the onset of hadronization, we may think of the effective degrees of freedom to be 
only (constituent) quarks with gluons essentially not being present as dynamical degrees of 
freedom. During hadronization, the effects of gluons will be primarily manifested in terms of 
confining color strings. In our picture the formation of hadrons will be with valence quarks 
which seems to be a reasonable approximation [9]. 

In the deconfined plasma, above T c , the color field lines are spread out, and quarks 
and gluons are quasi-free with the usual Debye screening of color charges. As the transition 
temperature is approached from above, confinement sets in and color field lines start squeez- 
ing in the form of color flux tubes connecting quarks with antiquarks forming mesons, or 
forming baryons, glueball etc. (We use here the language of color string for the simplicity 
of physical picture. Basic physics we use is simply that hadrons form by combination of 
nearby quarks/ antiquarks, as in the recombination models [9,10]). 

For simplicity, let us first consider a single quark and antiquark with the overall quark- 
antiquark system being in the color singlet state. Let us further assume a hypothetical 
situation when the quark and the antiquark do not have any kinetic energy as the system 
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approaches the transition stage. The color lines will squeeze in the form of a flux tube, 
connecting the quark with the antiquark thereby forming a meson. If the color field distri- 
bution about the line joining the quark and the antiquark was not very asymmetric then 
the resulting meson should most likely form in a state of zero orbital angular momentum. 
This is because the color field lines will just squeeze, forming the color flux tube, which will 
lead to the confining potential for forming the meson. This flux tube, for a symmetric con- 
figuration will not have any orbital angular momentum, but may have radial excitations, as 
it relaxes to the appropriate meson state in a given energy eigenstate. Even in a bag model 
like picture, initially static quark and antiquark will be attracted to each other, get bound 
inside the bag with appropriate kinetic energies inside the bag, developing the constituent 
masses. Again, for isolated pair, and symmetric initial color fields, it seems unlikely that an 
orbital excitation can result. 

Orbitally excited states for this case also can form due to following types of processes. 
As the quark-antiquark pair, connected by the flux tube (or inside the bag, in the bag 
picture), relaxes to an appropriate meson state, it can create some other lighter meson 
by fragmentation (or emit soft gluons if the initial quark and antiquark are not in a color 
singlet state, as in the Color Octet models [11]), resulting in the generation of orbital angular 
momentum for the remaining meson. Secondly, most often the initial color field configuration 
will be asymmetric (time dependent), leading to the formation of color flux tube (or bag) 
at the transition stage which is of deformed (time dependent) shape. Such a flux tube 
(bag), while relaxing to the meson state may generate orbital angular momentum. Also, 
orbital angular momentum can be generated due to the effects of other partons on this 
quark-antiquark pair. 

So far we have made the unrealistic assumption that the quark and the antiquark under 
consideration had no kinetic energies at the onset of hadronization. Let us now take into 
account the fact that at a non-zero temperature T, quarks and antiquarks will have thermal 
energies, with average kinetic energy per degree of freedom being about |T. For relativistic 
quarks one should talk about the total energy. However, at the onset of hadronization, 
quarks will develop constituent masses, so even for light quarks, the use of non-relativistic 
terminology may not be a serious problem, especially as our discussion is mostly qualitative 
in nature. Let us now consider the binding of such quarks and antiquarks into mesons at 
T = T c . The flux tube developing between the quark and the antiquark pair, now connects 
two point particles which have non-zero velocities. Clearly, in this case, for certain fraction 
(of order one) of cases, the quark and antiquark velocities will end up in the form of non- 
zero angular momentum of the resulting hadronic state. For example (using a semiclassical, 
non-relativistic picture), when their velocities are transverse to the flux tube connecting 
the pair, and are opposite to each other, then the resulting meson state will have non-zero 
orbital angular momentum. In the bag picture, with appropriate impact parameter, initial 
velocities will lead to orbital excitation as quarks get confined into hadrons. 

Note here that quarks will anyway develop appropriate kinetic energies as they fall into 
any bound state. Such kinetic energies will result from the color forces binding the quarks 
into hadronic states. As we mentioned above, this will happen even if quark/antiquark did 
not have any initial kinetic energies. However, if each quark and antiquark in the entire 
system has a given (average) definite kinetic energy, then in certain fraction of the cases 
when these velocities are appropriately oriented, orbital excitation should result. What 
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we are saying is that it is difficult to form L = hadrons when quarks/ antiquaries have 
large initial velocities, just as it is difficult to form orbitally excited hadrons with initially 
static quarks and antiquarks. L = meson can result only for special situations when the 
momenta of the quark and the antiquark are either parallel, or antiparallel with zero impact 
parameter. For generic situations, orbital excitation should result. (In the string picture, 
the string itself may have non-zero orbital angular momentum due to time dependence of the 
initial color field between the quark and the antiquark, as we mentioned above. However, it is 
not clear how, or whether, this string angular momentum may be related to thermal energies 
of partons, gluons in this case. Presumably, there will be no net orbital angular momentum 
left for the string if many gluons participate during the development of the color flux tube 
between the quark and the antiquark at the onset of confinement. Any other contributions 
to string angular momentum should be present even when the system temperature remains 
below T c , hence it should not affect the qualitative aspects of the signal we are proposing.) 

From the above intuitive reasoning, we conclude that thermal velocities of quarks and 
antiquarks at the hadronization stage can lead to the formation of orbitally excited mesons. 
In this picture, the thermal kinetic energies of the quark and the antiquark provide the 
required excitation energy for creating the orbital excitation. Now, in general the kinetic 
energy available will not be close to the required orbital excitation energy. For example, 
suppose the separation between the L = state and the orbitally excited state is of order 
800 MeV or more (as for the orbital excitations of n and rj mesons). Even the combined 
thermal kinetic energies of quark and the antiquark should not be much larger than about 
2-3T c ~ 500 MeV (with T c ~ 170 MeV). In such a situation the required orbital angular 
momentum has to come from other sources as well, as discussed above. 

However if one considers those specific mesons for which the orbital excitation energy 
roughly matches with the combined thermal kinetic energies of the quarks, then the required 
orbital excitation energy is right away available for creating the excited meson state. As 
every quark/ antiquark has an average kinetic energy of order T c at the stage when they get 
confined to hadrons, we conclude that formation of those orbitally excited mesons which 
have orbital excitation energies of order aT c should be greatly enhanced. Here the factor a 
is introduced to parameterize the part of kinetic energy of the quarks that gets converted 
into the orbital excitation energy of the meson. We are assuming here that a significant part 
of thermal kinetic energy of quarks survives as quarks develop constituent masses at the 
onset of confinement, even for light quarks. We again mention that quarks, when they get 
bound into hadrons, will have appropriate kinetic energies (e.g. in the bag model), which 
will result from confining color forces. What we are suggesting is that it is possible that a 
part (characterized by a) of initial thermal kinetic energies of quarks, when quark/ antiquark 
velocities are appropriately matched, may end up directly in the form of orbital excitation. 

With this, we can have a rough estimate of the range of a by first assuming that it is only 
one specific velocity component of, say, quark which will give the required orbital excitation 
energy, when combined with opposite velocity component of roughly similar magnitude of 
the partner antiquark. Each velocity component will have thermal kinetic energy of about 
|T C at the transition stage (again, using a non-relativistic picture with quarks which develop 
constituent masses). Thus the total thermal energy of the quark-antiquark system which 
can be used to give orbital excitation energy will be about T c . Other velocity components 
of quark and the antiquark could combine to give higher orbital excitation energy, or could 
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contribute to radial excitation energy. For example, average kinetic energy of the quark 
for the degree of freedom transverse to the color string will be about T c . Combined kinetic 
energy of the quark and antiquark for suitable configuration (with opposite velocities for 
quark and antiquark) can then become 2T C which can directly translate into the orbital 
excitation energy. From the above qualitative discussion we expect that the average value 
of a could be close to 1-2, and in any case should not be larger than 3. The value a = 3 
corresponds to the case when entire kinetic energy |T of the quark and of the antiquark 
ends up in the form of the orbital excitation energy, (of course, these are kinetic energies 
on the average only.) For the case of baryons, the estimates of a should be appropriately 
modified by accounting for the thermal kinetic energy of each of the three quarks. 



III. EXCITED HADRONS IN RELATIVISTIC HEAVY ION COLLISIONS 

Let us now consider the general situation of a relativistic heavy-ion collision. Let us 
first consider the situation when the collision process creates a system which thermalizes 
but does not reach upto the transition temperature T c so the QGP state is never created. 
In such a situation the formation of hadrons can be modeled using string fragmentation as 
the QCD strings break creating secondary hadrons [12]. In this situation there appears no 
reason to relate the thermal energies (of hadrons) to the orbital excitation energies, at the 
stage of hadron formation. The production of orbital excitations will mostly result from 
the transverse momentum created for the quark- antiquark pairs from the breaking of QCD 
string in the fragmentation model [12]. There will be a distribution of P T of quarks and 
antiquarks, with no specific value of Pt preferred. In particular, there is no reason to relate 
the transverse momentum of quarks/antiquarks at the fragmentation stage to the thermal 
energies in the resulting hadron system. Thus there is no reason to expect enhancement 
in the multiplicities of any specific orbitally excited states at some particular value of the 
temperature of the final hadronic system. 

Now consider the situation when the temperature of the system reaches above T c and 
deconfined QGP state forms. As the QGP expands and cools, its temperature drops to T c 
and partons combine to form hadrons. As we argued above, now specific orbitally excited 
states, with orbital excitation energies close to aT c will form with enhanced rate. We should 
note that the formation of hadrons in a QGP state which is cooling through T c and the 
formation of hadrons in collisions where temperature never reaches upto T c may be of very 
different nature [9,10]. Therefore one may not be able to compare directly the multiplicities 
of orbitally excited states between the two cases. However, we can consider the following 
ratio, 

Here N(M*) and N(M) are the total multiplicities of orbitally excited meson and the 
corresponding L = state meson respectively. For such a ratio, the effects of any possible 
intrinsic difference between the two types of hadronization may be minimized. Also, to make 
sure that hadrons which form from hadronization of QGP do not get mixed with hadrons 
forming in the regions which do not undergo deconfining transition one should confine to 
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near central rapidity regions. As we have mentioned above, restriction to hadrons with not 
very high should ensure that only hadrons forming via a recombination type process are 
being included and not the ones coming from the fragmentation process. 

As for other orbitally excited states, for which the orbital excitation energy is very 
different from the expected range of aT C) for them also the ratio R* defined above may 
be somewhat enhanced due to thermal velocities of quarks and antiquarks. However, for 
them the dominant contribution to the orbital excitation energy must come from sources 
other than the thermal energy. Those sources will lead to contributions which may become 
difficult to distinguish from various sources for the case of hadronization when temperature 
remains below T c . Hence, in those cases the ratio R* may not show any significant change 
as the maximum temperature of the system crosses T c . 

It is not entirely clear what happens for the states which have excitation energies well 
below T c . Though, in the particle data table there appear to be no such states, still, con- 
ceptually one should discuss this case. If thermal kinetic energy of quarks and antiquarks 
is very large then during the formation of hadrons the remaining kinetic energy has to be 
either converted to part of the energy of QCD string, or dissipated away in other modes, 
e.g. by emission of lighter mesons. 

The effects we have discussed above are valid for baryons as well, as the basic physics 
of orbital excitation coming from the transverse components of thermal velocities of quarks 
remains the same. It is also possible that the effect we are proposing here may be present 
even for radially excited states of mesons. However, radial excitations may not be uncommon 
even when quarks do not have significant thermal velocities. For example, when a quark 
and an antiquark become joined by the QCD string at the stage of hadronization, then 
while relaxing to an appropriate meson state, radial excitations can result. It may thus 
become difficult to disentangle the effect of thermal kinetic energies on radial excitations. 
We therefore focus on the orbital excitations of hadrons as the underlying physical picture 
is simple to analyze. 

The physical picture which we are using here for the enhanced formation of specific 
orbitally excited states is similar to the formation of resonances where the cross-section 
shows a peak when the center of mass energy equals the energy of the intermediate bound 
state. The difference here is that we are here focusing primarily on the initial kinetic 
energy of the quark and antiquark and essentially arguing that the cross-section should show 
enhancement when the initial kinetic energy is of the same order as the orbital excitation 
energy of the meson state. Clearly this special treatment of initial thermal kinetic energies 
of quarks is a very strong assumption, even stronger than the non-relativistic, semiclassical 
picture underlying the entire discussion. For heavy quarks the non-relativistic treatment is 
reasonably accurate [13]. Even for light quarks, as the quarks start binding into hadrons 
the constituent quark mass will start playing the dominant role and the non-relativistic 
treatment may still capture essential physics. (Velocities of light quarks are about 0.7-0.8 in 
heavy- light bound states [14].) What remains a gross approximation is the special status of 
the initial (thermal) kinetic energies of quarks/antiquarks. As we discussed above, quarks 
and antiquarks, when they get bound into hadrons will develop kinetic energies (e.g. in the 
bag model), which will result from confining color forces. What we have focused on is that 
the system undergoing hadronization consists of quarks and antiquarks (as we mentioned 
earlier, we consider only on the valence quarks/antiquarks in hadron formation), all of which 
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have roughly same kinetic energies (of order T c ). When such rapidly moving quarks and 
antiquarks get captured into hadrons, then it seems very natural to assume that it will be 
relatively difficult to form L = hadrons and that a part (characterized by a) of initial 
thermal kinetic energies of quarks may end up directly in the form of orbital excitation. 

As confinement sets in, and quarks pick up constituent masses, the kinetic energies of 
quarks will contribute to resulting potential energy of the system, e.g. to the energy of 
flux tube which connects a quark and an antiquark to form a meson. Kinetic energies of 
quarks will also get significantly affected as they develop constituent masses due to effects 
of confinement. Such effects should tend to lower the value of a, and clearly all such effects 
have to be properly accounted for. For example, the spin-orbit interaction energy also may 
play an important role in determining the value of a. In the absence of a detailed picture 
of confinement, we just parameterize such effects in terms of a. 

Even with all these ambiguities and crudeness of our arguments, we suggest that there 
remains the possibility that the thermal velocities of quarks and antiquarks may have a siz- 
able contribution to the orbital excitation of hadrons (and possibly also for radial excitations 
of hadrons). If this is true then the relative ratios of orbitally excited states R* (as discussed 
above) carry a signature of the transition temperature T c at which all the quarks/antiquarks 
become bound into hadrons. 

The specific signal we propose is the following. One should plot the ratio R* for various 
hadronic states as a function of the center of mass energy of the collision, for, say, given 
colliding nuclei (which determines the initial energy density, and hence the temperature of 
the resulting system [15]). Alternatively, one can also plot R* as a function of number of 
participants N part /centrality for a given center of mass energy. For generic hadronic states, 
R* should vary in a certain smooth manner as a function of y/s (or, of centrality). For 
example, if we assume that abundances of hadrons are determined by the thermal equilibrium 
distribution, as seems to be the case for most hadrons [8,16], then the dominant behavior of 
the ratio R* should be given by, R*h erma i ~ exp[— (E* — E)/T] (with appropriate chemical 
freezeout temperature). Here E* and E are the energies of the orbitally excited hadron and 
the L = state hadron respectively. We are not writing here the detailed expressions for 
the multiplicities (including possible factors of the appropriate chemical potentials), as in 
the thermal models [8]. Also, we do not discuss about the contributions of decays of heavier 
states to various multiplicities at freezeout. This is because such factors will be present for 
the ratio R* even when the system temperature always remains below T c , and hence do not 
affect the nature of the signal we propose. For our purpose it suffices that R* is some smooth 
function of y/s. 

It is important to clarify here that we are mentioning this case of thermal abundances 
only as an example. We are not arguing for any specific shape of the plot of R* vs. y/s (or, 
Np a rt). For example, it is very likely that resonances do not chemically equilibrate, so the 
assumption of thermal distribution for them is inaccurate. Further, the final multiplicities 
will be severely modified by factors such as final state interactions and loss of signals in the 
medium, which will be expected to strongly suppress resonances as a function of y/s (or 
centrality), (also by enhancement factors such as resonance regeneration). The overall curve 
of R* may very well be decreasing as a function of y/s or N part (as for A*(1520)/A, see ref. 
[6,7]). 

Now, at some critical value of y/s at which the temperature of the system reaches 
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above T c , the above described enhancement of orbitally excited states should become oper- 
ative due to the formation of QGP state. This should result in an abrupt enhancement in 
the value of R* for those hadrons for which the orbital excitation energy is close to aT at 
T = T c , i.e. 

E* -E~ aT c . (2) 

Further, this enhancement should prevail at all higher values of y/s. Again, note that 
if the overall curve of R* is decreasing as a function of y/s then our argument implies that 
beyond a critical value y/s~ c (or centrality), there should be a separate segment of the plot of 
R* with a jump enhancement at y/s^. Even for other hadrons for which the orbital excitation 
energy is not too close to aT c , there may be some enhancement in the corresponding R* 
beyond y/s~ C) with the enhancement decreasing as E* — E deviates more and more from the 
value of aT c . One can also consider the ratio of multiplicity of an orbitally excited hadron 
state, say with energy E*, to L = state of some different hadron which has the same (or 
almost the same) energy. Such a ratio, from thermal models, should not strongly depend on 
temperature (assuming other factors like chemical potential to be same for the two states), 
and hence on y/s, apart from the effects of final state interaction in the hadronic medium. 
Our model will then predict a sharp change in the value of this ratio at some specific value 
yfs~ c of y/s. This may have an advantage that departure from a (possibly) roughly constant 
value of the ratio at may be easier to identify. On the other hand, using states of 
the same hadron is advantageous in minimizing the effects of details of the parton system 
evolution etc. 

Note that this requires that the enhanced multiplicities of these specific orbitally ex- 
cited states should depart from the thermal distribution. This appears to be the case for 
resonances. In any case this possibility does not seem unreasonable to us as, even with 
the general success of statistical models, it is clearly a very strong assumption that strictly 
all hadronic multiplicities should follow equilibrium distributions. For this to happen, the 
hadronic system should last for long enough time (compared to all relevant hadronic inter- 
action rates) before the chemical freezeout happens. This does not seem to be the case as 
chemical freezeout temperatures do not seem to be much different from the expected value 
of the transition temperature [8]. In fact, it has been argued in ref. [8] that, in view of 
the short time scale between the hadronization and the chemical freezeout, hadronic multi- 
plicities are most likely established during hadronization. As the detailed dynamics of this 
process is not understood, is not clear to us how (or whether) this can lead to a suppression 
of the signal we are proposing. It is possible that if there is a large enough contribution to 
the multiplicities of certain specific orbitally excited states, as we are proposing, then the 
resulting jump in R* may survive. 

We do not know the value of a so one cannot say exactly where this enhancement should 
occur. Still at some value of T, and hence of y/s (or centrality), the enhancement in the 
value of R* as described above should occur. Such an abrupt enhancement will signal the 
formation of QGP, and at the same time, the value of T (extracted from the value of y/s) 
[15] will give the value of the transition temperature T c if we can control the estimate of a. 
As we argued above we expect a to be roughly in the range of 1-3, with more probable values 
being close to 1-2. Thus the value of aT c can be conservatively estimated to lie somewhere 
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in the range of 150 - 500 MeV, with a value of about 150-350 MeV being a more probable 
estimate. 

We emphasize here three important aspects of the signal we have proposed. First, a 
jump in the values of R* at a common value y/sl of y/s (or of centrality), for a set of 
hadrons, will be a concrete signal that the system has gone through a transient stage of 
QGP phase. Secondly, from the value of y/s~ c one can extract the value of T c depending on 
the value of a, as mentioned above. This determination of T c , in some sense, will be similar 
in character to the determination of T c through the value of y/s at which the proposed J / ip 
suppression [2] takes place due to QGP formation. (Though, note that J/ip suppression 
happens at a significantly higher temperature than T c , above about 2T C , see ref. [17]. Also, 
J/ip suppression is not expected to be sharp at some definite y/s, it gets smeared because 
of final state interactions. As there appears no reason to expect such strong smearing 
in our model, the jump in R* may be sharper.) Apart from the uncertainties in a, this 
determination of T c will also depend on the details assumed for the entire collision process 
which relate the center of mass energy of the collision to the energy density of the plasma, 
and hence temperature by assuming thermal equilibrium [15]. Various parameters, such 
as the initial thermalization time, crucially affect these estimates of the temperature. The 
third feature of our proposed signal is important from this point of view, as here one will 
have another completely independent determination of the approximate value of T c . This 
arises from the fact that the specific orbitally excited states of hadrons which will show the 
jump in corresponding values of R* are precisely those which have orbital excitation energies 
of order aT c . For example, the candidate states we have proposed (see below) have been 
chosen because the expected value of T c is about 170 MeV. (From this point of view, one 
should plot R* vs. y/s for all available partner hadronic states, to find out which ones, if 
any, show a jump in R*.) Of course, as the value of a is uncertain at least by a factor of 
2, one will only be able to infer the value of T c within that factor by looking at the orbital 
excitation energies of states which show jump in R*. Still, it is important to appreciate that 
this determination of T c is completely independent of the determination of T c through y/s~ c 
(or specific centrality) at which the jump in R* is expected to occur, and is uncertain only 
by factor a. (Of course, for better estimates of a one needs a proper understanding of the 
hadronization process itself.) 

We mention another important effect which naturally happens in our model. In our 
picture, the thermal kinetic energies of quarks get used up in generating the orbital excitation 
energies for certain specific hadrons (as given above). Thus, for these hadrons, immediately 
after formation, there should be little thermal kinetic energy left for the center of mass 
motion of the hadron. For other hadrons for which the thermal kinetic energies of the 
constituent quarks does not result in orbital excitation, net kinetic energy of hadron should 
be significant. Thus, one expects that orbitally excited states should have low Pt, while the 
L = state hadrons should have Pt as consistent with the thermal kinetic energies involved. 
We will use this fact in the estimation of R* using a recombination model in the next section. 
As the kinetic energy of the order of aT c gets used up in creating the orbital excitation, we 
expect that, on the average, Pt of the excited states should be determined using kinetic 
energies which are lower by about aT c compared to other states. Of course, this difference 
exists only at the stage of hadronization. If thermal equilibrium of hadrons prevails for 
any significant time after hadronization, all hadrons will acquire thermal kinetic energies. 
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However, for certain collisions it is possible that hadronic thermal freezeout happens almost 
immediately after hadronization [20]. (For example, with large value of \/s to allow for 
QGP formation, but with small nuclei, so that transverse expansion becomes dominant at 
an early stage.) If that happens then orbitally excited states (or, in general even the radially 
excited states) should have significantly low P T compared to other hadrons. (In principle, 
this effect may be operative for even those orbitally excited states which have significantly 
larger excitation energies compared to aT c , as part of the orbital excitation energy may 
still be coming from the thermal kinetic energies of quarks.) These arguments have the 
important implication that for such collisions (where thermal freezeout happens close to 
hadronization), the jump in R* which we have proposed to occur at some critical value of 
y/s, should be much more prominent if one restricts to hadrons with low Pt- 



IV. ESTIMATE OF R* IN RECOMBINATION MODEL 

Some estimate of the ratio R* can be obtained within the framework of recombination 
model discussed in ref. [9]. For rough estimates, we first review the non-relativistic estimate 
of the momentum distribution of the ground state meson M from ref. [9]. Though we 
have quoted data for A baryon from RHIC ( [6,7]), we will not attempt to calculate R* for 
baryons in the recombination model to keep the basic physical arguments simple. After all, 
recombination models do not capture the physics of development of confining forces (string) 
between quarks and antiquarks. Therefore, the calculations provided in this section still 
have to be supplemented with the basic physics of our model, namely, certain situations will 
be argued to lead to L = meson while others leading to L — 1 meson. 

Spatial wave functions for a (quasi) free quark and antiquark (with momenta pi and P2), 
and for a meson (with momentum P) formed by this quark- antiquark system are written as, 

< x\Q, pi, p 2 >~ e i(pi xl+P2 X2) , 

<x|M,P>~e lP - R M (y). (3) 

Here, R = (xi + x 2 )/2 is the center of mass coordinate and y = xi — x 2 is the rela- 
tive coordinate. 0m (y) is the normalized meson wave function. The relative momentum 
conjugate to y is q = (pi — p 2 )/2, and the center of mass momentum is P = pi + p 2 . 

Meson momentum distribution is then evaluated using the overlap < Q,pi,p 2 |M, P > 
(see, ref. [9] for details). One gets, 

Here, tu(q) gives the parton phase space distribution, 0m (q) is the meson wave function 
in momentum space, and A accounts for constant factors like degeneracy and volume etc. 

Let us consider the situation when the momenta p x and p 2 of the quark and the antiquark 
are large (compared to the width of 0m(q), which will be expected to be of order Aq C d)- 
We then consider two situations, 

(1) pi and p 2 are almost parallel, as discussed in ref. [9]. Then |q| << |P|. As we 
discussed, this situation will correspond to the formation of L = meson, with the meson 
center of mass momentum P being large. 



11 



(2) pi and P2 are almost opposite to each other. Then P is small and |q| >> |P|. This 
is the opposite limit of the recombination case discussed in ref. [9], where this situation is 
not considered as leading to formation of any meson. For example, in ref. [9], it is argued 
that hadrons with large P can result from combining quarks with smaller momenta. We 
take this to hold true for L = meson, as mentioned in the case (1) above. However, for 
L = 1 hadrons we take the opposite situation where quarks with large momenta combine to 
give a hadron with small P. In ref. [9], higher L states are incorporated by using degeneracy 
factors, as in the statistical models. In contrast, in our model, almost opposite momenta of 
quark/ ant iquark result in an orbitally excited meson, (with appropriate impact parameter 
between the quark and the antiquark). Note that with long range confining forces operating 
between quarks/ ant iquarks during hadronization, combining even fast moving quark and 
antiquark (with opposite velocities) to form a meson should not be a problem. We can 
then incorporate this basic difference to modify the approach of ref. [9] for discussing L = 1 
meson. 

We first repeat the main steps of case (1) from ref. [9], and then modify those steps 
for case (2). For case(l), the meson multiplicity is calculated ( [9]) by using the Taylor 
expansion of u>(^ ± q) about y (as |q| << |P|). This gives, 

W (| + qM|-q)= W (|) 2 + 0(q 2 ). (5) 

With this, the lowest order contribution to the meson multiplicity can be obtained ( [9]) 
just by using the normalization of the wave function of the meson. One gets, 

dN M P 2 __p! 

-— = Aw(-y = Ae ^t. (6) 

Here we have used the non-relativistic distribution for the quarks with mass m (for 
simplicity we use same flavor for quark and the antiquark), 

P 2 

W(P) ~ e 2mT. (7) 

Important thing to note in Eq.(6) is that to this order the expression for the meson 
distribution does not depend on the shape of the meson wave function, which, for the 
ground state meson can be taken as (in momentum space) 

<Mq)~e^ 2 / 2A M. (8) 

Am is the width of the wave function, and should be of the order of Aqcd- 
Now we consider case(2) and discuss the formation of L — 1 meson. Since for this case 
we consider |P| << |q|, we Taylor expand u>(y ± q) about ±q. We get, 

™(q + 2>(-q + 2") = ™(q>(-q) + o(P). (9) 

Neglecting O(P) terms, we get, 

w =A 7(^ w(qM - q)| ^ (q)|2 - (10) 
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A' differs from A in Eq.(4) in the degeneracy factor. It is clear that in this limit 
(|P| << |q|), the form of the wave function will be important in determining the mo- 
mentum distribution of the meson M*. Since in this case we are considering formation of 
L — 1 meson, we use the following form of the wave function for meson M* ( [18]), 

0m* (q) ~ q x e (11) 

Note that we use the same width as for L = meson. Here, q x is the x component 
of q. Using Eq.(7) for iu(±q), we get, 

dN M * f d 3 q _qi 2 

A / TTT-^e »»Tg2 e m . (12) 



Let us define, 



Then we can write, 



d 3 P 7 (2tt) 



q' = q(l + ^k)V2 = qC7 . (13) 
ml 



dN M * A' f d 3 q' /2 -ft- , . 



d 3 P (2tt) j 

The integral is unity due to the normalization of the wavefunction. We thus get, 

dN M » _ A' 
d 3 P ~ C E ' 



(15) 



Using this equation and Eq.(6) we get the following expression for the (momentum 
dependent) ratio R*. 

p 2 

R * P = (l + AVrnT)^' (16) 

with T here being equal to T c . Here, k gives the ratio of the two degeneracy factors. 
This expression shows a rapid increase for large value of [L = 0) meson momentum. One 
may consider this increase as coming from the trivial factor of e _p / 2mT for the momentum 
distribution of the L = meson, as the L = 1 meson is taken here to have small momenta. 
However, it is important to note here that both mesons (L — and L — 1) are considered 
here as forming from the same quark and antiquark having large momenta pi and p2- There 
is a suppression factor for both, quark, and antiquark, of the form e~ Pi / 2m , i = 1,2. For 
L = meson, this leads to large momenta P for the meson, leading to the suppression 
factor e~ p l 2m . On the other hand, for L = 1 meson, the same suppression factors of the 
quark and the antiquark are absorbed in the wave function of the meson. This softens the 
suppression factor from an exponential suppression to a multiplicative factor as in Eq.(16). 
Note that this expression should eventually breakdown at sufficiently high P. This is because 
our assumption that a quark and an antiquark of opposite momenta leads to an (orbitally 
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excited) bound state will not hold true for very large values of quark/ ant iquark momenta. 
Even for L = meson, as discussed in ref. [9], recombination model is not expected to 
be applicable for very large values of P (beyond a couple of GeV), as the physics may be 
dominated by fragmentation in this regime. 

We should clarify that we are not claiming here a derivation showing the enhancement of 
R* as we have argued in previous sections. For example, one could have very well considered 
the L = meson wave function for case (2), instead of Eq.(ll), with similar result as in 
Eq.(16) with the exponent 5/2 changing to 3/2. (Or, one could have considered other 
excited states of the meson.) The above estimate of R* P has to be supplemented with our 
basic picture that for a quark and an antiquark with large momenta pi and p 2 , L = 
meson is likely to result when the two momenta are almost parallel. On the other hand, 
if the two momenta are almost opposite to each other, then L = 1 meson (or, possibly a 
radially excited one which we do not discuss) is likely to form because the available thermal 
kinetic energy may be enough to lead to L — 1 excitation (at least for the candidate states 
we suggest in the next section). The estimate of R* provided here is, therefore, to be taken 
as indicative of the extra contribution to the multiplicity of the orbitally excited meson M* 
when QGP formation occurs. 



Using the particle data table [19], there are several candidate hadrons which have or- 
bitally excited states with excitation energies of order 200-400 MeV above the corresponding 
L = states. We give several such candidates in the following. Other candidate states can 
be found from ref. [19] with the basic criterion being the relevant range for the orbital exci- 
tation energy. Below, each equation gives the hadronic state, along with its partner orbitally 
excited state. 



V. CANDIDATE HADRONIC STATES FOR THE SIGNAL 



Mesons : 




l - ) m = 2112 MeV 
+ ) m = 2317 MeV 
Am = 205 MeV 



(17) 



D*°(J P 

d;°(j p 



r) m = 2007 MeV 
2+) m = 2459 MeV 
Am = 452 MeV 



(18) 



J/^(1S)(J P = r)m = 3097 MeV 
Xco(lP)(J P = + )m = 3415 MeV 
Am = 318 MeV 



(19) 



Baryons : 
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A C {J 


= 2 } 


m = ZZoo MeV 




"2 } 




k+{J P 


= 1") 
2 


m = 2625 MeV 




Am 


= 309, 340 MeV 


A ( T p 


= 2 } 


m = lllo MeV 


ill J 


"2 } 


™ — 1 jnfi MpV 
/7 6 — l^JzUU IVlt; V 


A(J P -- 


9 


m = 1519 MeV 




Am = 


= 290,403 MeV 


~ + ( i p 


"2 } 


m = z4oo MeV 


~UJ p 


3 + 

= - ) 

2 ; 


m = 2647 MeV 






m = 2790 MeV 




Am 


= 324, 143 MeV 



(21) 



(22) 





-r> 


m = 2472 


MeV 


~Uj p 


3 + 

= - ) 

2 ; 


m = 2645 


MeV 


~UJ p 




m = 2790 


MeV 




Am 


= 318,145 


MeV 



We mention that for some of these states J as quoted above are not confirmed, see 
ref. [19]. The energy differences for the partner states as given above are reasonably close 
to the expected value of aT c (with T c ~ 170 MeV), and hence should be good candidates 
for evaluation of the ratio R* as described above. For example, £>^(2317)(J P = 0+) and 
D* S ± (2112)(J P = 1") have a mass difference of 205 MeV and hence one may expect that an 
important fraction of the orbital excitation energy in this case may come from the thermal 
kinetic energies of quark and antiquark (even with lower estimates of a as mentioned above). 
In principle one could allow the spin of the L = state and the corresponding orbitally 
excited state to be different. In view of hyperfine contributions one may like to have the same 
spin for the partner states for the evaluation of the ratio R*. However, such contributions 
will be present for the formation of these states even below T c , and for hadrons forming at T c 
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in our picture it may just affect the vaiue of a without affecting the qualitative aspect of the 
signal (i.e., a jump in the value of R* at y/fQ. We have not addressed here the important 
issue of experimental observations of all these states mentioned above, and for many of these 
it may be difficult to have data in relativistic heavy-ion collisions. 

Recently, measurements of A*(1520)/A have become available from RHIC as a function 
of N part [6,7]. Note that this set is one of the candidates we have listed above (Eq.(21), with 
A* listed at 1519 MeV). The multiplicity ratio R* for A*/A shows a strong decrease as a 
function of N part [6,7]. As discussed above, this overall decrease is of no consequence to our 
model. However, the value of the multiplicity ratio for very large N part (beyond 350) may 
be consistent with a slight increase, though completely within the error bar, to conclude 
anything definite. We again stress, what is important is that if there is any increase in 
R* beyond some N pari then it will be very hard to explain it with existing models where 
a rapid decrease in the multiplicities of resonances arising from final state interaction in 
the hadronic medium is naturally expected. Our model will naturally account for any such 
increase. In fact, the (possible) increase in R* seems to occur beyond sufficiently large 
value of N par t ~ 350 which may be consistent with QGP formation. We mention that this 
(possible) rise in R* for A*/A is more distinct in ref. [6] where R* is plotted as a function of 
N pari . In ref. [7], the plot is w.r.t. dN ch /drj, and the rise does not appear clear (which may 
be due to different binning). One has to wait for more refined data with a more detailed plot 
of R* vs. N part to see if actually there is a discontinuity in the plot of R* for A*(1520)/A 
beyond some value of N part . 

It is important to realize here that the abundances of the hadrons may be dependent on 
QGP formation from other physical considerations also. For example, the suppression of J/ip 
and Xci due to Debye screening of color in the QGP state is well discussed in the literature. 
Further, due to larger radius of Xc t , the suppression from Debye screening is larger for \ Ci 
[17]. Even in the recombination model, one would expect that forming orbitally excited 
states may be harder due to their larger sizes as the medium is very dense. Clearly such 
effects will give rise to opposite behavior for the ratio R* as compared to the effect of thermal 
kinetic energy of quarks as we have discussed. Final dependence of R* on y/s will result from 
a combination of the two effects. (Note that similar effects also occur from hadronic final 
state interaction. However, such hadronic state interaction effects are of no consequence to 
our model which only predicts a jump in R* beyond some y/s, and not the overall shape 
of the curve of R*. Thus, decrease of A*/A vs. N part may be naturally expected from such 
final state interactions, but any possible rapid rise beyond some value of N part cannot be 
accommodated in such models.) Though there are detailed estimates [17] of suppression 
of J ftp and x^, it is not possible to say which effect will dominate due to lack of definite 
quantitative estimates for the effect we are proposing. In any case, for hadrons with light 
quarks, the Debye screening is very effective due to larger sizes of these hadrons (compared 
to the Debye length), so one only needs to consider their formation from quark/ ant iquark 
recombinations at T c . Thus, for these hadrons, the predicted jump in R* should not be 
adversely affected by these other considerations. 
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VI. CONCLUSIONS 



In conclusion we have suggested a novel signature for the detection of the transient QGP 
state formation in relativistic heavy-ion collisions. The basic nature of the signal is almost 
model independent in the sense that it does not depend on the nature or even the existence 
of phase transition between the confined and the deconfined state of QCD. Further, details 
of plasma evolution etc. also may not be of much importance. This is because R* depends on 
the ratios of multiplicities of states of a given hadron and presumably the effects of plasma 
evolution etc. will get minimized in taking the ratio. 

We emphasize that though the arguments we have presented are very qualitative and 
crude, the basic physics underlying these arguments is simple. Main argument we have used 
is that quarks/antiquarks have thermal velocities at the stage when they all get captured 
to form bound states at the hadronization stage, i.e. when the temperature of the system 
is equal to T c . Due to this thermal velocity, when QCD string starts joining, say, a quark 
with antiquark to form a meson, it is easy to generate orbital motion of quarks in the bound 
states. Essential part of the argument being that it is difficult to form L = hadrons 
when quarks/antiquarks have large initial velocities, in the same way as it is difficult to 
form orbitally excited hadrons if quarks/antiquarks did not have any initial velocities. We 
then further argue that for those hadrons for which the orbital excitation energy is a certain 
fraction (characterized by a) of the thermal kinetic energy of quark/ antiquark system at 
T = T c , the formation of such orbitally excited bound states will be enhanced. Using these 
arguments we suggest that the ratio of multiplicities of orbitally excited state to the L = 
state of such specific hadrons should show an abrupt increase at a critical value of \/s which 
corresponds to the formation of an intermediate, transient, QGP state. As we mentioned 
above, alternatively, one can also plot R* as a function of centrality for a given center of mass 
energy, with R* expected to show a jump above certain value of the number of participants. 

The signal we have proposed is simple and just requires plotting of the respective ratios 
R* (as defined by Eq.(l)) for various hadrons with respect to ^/s. One important feature of 
this signal is that it predicts a jump in the value of R* , at a common value of (for given 
colliding nuclei), for a whole class of hadronic states which have orbital excitation energy 
splitting of similar magnitude. Most likely candidate hadrons which have orbital excitation 
energies roughly in the required range are charmed and strange hadrons. For example, 
L>^(2317)(J P = 0+) and D* S ± (2112)(J P = 1") have a mass difference of 205 MeV. Thus, 
a plot of the ratio of total multiplicities of L>^(2317)(J P = 0+) and D* S ± (2112)(J P = 1") 
with increasing ^/s may contain the information of QGP formation as discussed above. We 
have given a list of several such candidate hadronic states. We note that the recent data [6,7] 
on the multiplicity ratio A*(1520)/A from RHIC as a function of N part may be consistent 
with our model, with a possible hint of a slight increase beyond N pari ~ 350. One needs to 
have a more refined data, with more detailed plot for intermediate values of N part to make 
any definite conclusions. We have not addressed the issue of experimental observations for 
various resonances listed in the previous section. It may not be easy to collect enough 
statistics for many of these candidate states [21,3]. As we mentioned above, the effect we 
have proposed may be present even for radially excited states of hadrons. However, for radial 
excitations the effect may not be prominent. We have primarily discussed orbital excitations 
of hadrons as the basic physics of our arguments is simpler to model for this case. A better 



17 



understanding of the underlying physics of quark/ ant iquark recombination into hadrons is 
needed before one can make any definitive quantitative predictions regarding this signal. 
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